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Метод дискриминантного анализа моделей 
пласта-коллектора по результатам 
гидродинамических исследований скважин 


В статье описывается новый количественный метод для дискриминантного анализа моделей, который был 
назван методом последовательного прогнозирования вероятностей. Данный мегод базируется на 
байесовском выводе. Метод является прямым продолжением способов, использующих концепцию 
доверительных интервалов, но в то же время лишен их основных недостатков. 


Метод последовательного прогнозирования 
вероятностей 


Идея метода последовательного прогнозирования вероятностей (МИПВ) состоит в 
том, что корректная модель должна более точно предсказывать изменения давления, чем 
другие модели. С количественной точки зрения вероятность реального изменения 
давления в будущие моменты времени для корректной модели должна быть выше, чем 
прогнозируемая вероятность для остальных моделей. Прогнозное распределение вероят- 
ностей событий, связанных с изменениями давления, может быть рассчитано при извест- 
ных распределениях вероятностей параметров. 

С другими методами, схожими с предлагаемым подходом, можно ознакомиться в [1]. 

Базовый алгоритм МИПВ для проведения дискриминантного анализа моделей 
пластов на основании гидродинамических исследований скважин (ГДИС) имеет сле- 
дующий вид. 

1. Выбрать несколько моделей-кандидатов пласта, которые так или иначе 
согласуются с данными ГДИС и другой доступной информацией. 

2. Использовать первые несколько результатов замера давления для оценки 
параметров пласта и предсказать распределение вероятностей для давления в 
следующий момент времени применительно к каждой модели пласта. 

3. Вычислить вероятность путем подстановки реального значения давления в 
формулу для прогнозируемого распределения вероятностей и обновить значение 
суммарной вероятности, умножив ее на только что рассчитанную вероятность, 
применительно к каждой модели пласта. 

4. Повторять шаги до тех пор, пока не будет достигнуто различие в суммарных 
вероятностях для каждой модели пласта. 

5. Дискриминантный анализ моделей-кандидатов пласта осуществляется на основании 
вычисленных суммарных вероятностей. 
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Плотность вероятности у, Плотность вероятности у, 
на основе Модели 2 на основе Модели 1 


Прогнозное рас- 
пределение веро- 
ятностей у, на 
основе Модели 2 а. 


Прогнозное рас- 
пределение веро- 
ятностей у. на 
основе Модели 1 


——— 
Текущие данные Будущие данные 


х 
Рисунок 1 — Схематическое представление метода 
прогнозирования вероятностей для двух моделей 


Прогнозное распределение вероятностей для каждой модели определяется в соот- 
ветствии со следующей процедурой в рамках байесовского вывода. При этом параметры 
и считаются случайными величинами с некоторыми распределениями вероятностей. 

Предположим, что нам известна дисперсия с” этих распределений. Тогда неоп- 
ределенность, связанная с моделью, описывается следующим образом: 

не ен 
РиоБ(иу,,.... у,)= те оч ре ( и) н("-*) (1) 


(2по’) 


Параметры и образуют многомерное нормальное распределение вокруг Я с 


матрицей ковариаций о’Н"' при известных и независимых наблюдений результатов 
измерения давления. и — это вектор размерности т, а значит, их распределение веро- 


ятностей также имеет размерность т. Средние (ожидаемые) значения параметров и 
являются оценками и по методу наименьших квадратов (МНК) на основе и незави- 
симых наблюдений. Матрица ковариаций о”Н`' — это тхт матрица. 

С использованием оценок параметров и и и наблюдений реальное значение 


давления у’, в точке х„, выражается путем разложения в ряд Тейлора с сохране- 


п-+1 


нием членов вплоть до первого порядка: 


Ун=Р(и, 1) = Е(и, ха) +8" (ии), (1) 
где 
и 
вы Е (и, х,.:) о ОЕ (и, х.1) | 6) 
00, 00, ие 


Градиент функции модели © вычисляется в точке х,., с помощью оценочных 


п-+1 
значений параметров на основе первых и наблюдений. 
Выбрав обозначение у. =Ё (м. а ‚ уравнение (1) можно представить в виде 


Уи-у,и= 8" (и-и). (3) 
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й С 
В соответствии с (1), У, нормально распределено относительно у,., с прог- 
нозной дисперсией 


2.2. 9-1 —_ Т-! 2 
= ОН -5=8Н 5.0, (4) 
и условной вероятностью 
* ыы * ИЕ 1 1 * — г 
РОБ (у, У, ) = РОБ о) про оао: н) 5“ № 


Уравнение (5) характеризует неопределенность, связанную с моделью, которое 
показывает, насколько она адекватна, и полностью определяется средним (ожидаемым) 


^ 


> > 2 
значением у,,., и прогнозной дисперсией с’. 
Неопределенность, связанная с данными, описывается следующим образом. 
Предположим, что наблюдаемое изменение давления у, нормально распределено 
относительно реального изменения давления у”, с дисперсией о’: 


това) ее Бу. %„} (6) 


Точно такое же допущение делалось при выводе (1). Формула (6) выражает не- 
определенность, связанную с данными, которая характеризует то, как хорошо модель 


соответствует данным. Поскольку точное значение у’. неизвестно, связь между у, И 


у„.: Устанавливается путем исключения при интегрировании по у’. : 


РБ, н|Уа| = Рио (у. у,}= 


= [РБ (ууу) Рио (у У) = (7) 


ар (рана) 
2т (с? +) 2(5* +5: 


В результате наблюдаемое изменение давления у, нормально распределено относи- 


тельно прогнозируемого (ожидаемого) изменения давления у„., с общей прогнозной дис- 
> - 2 

персией с” + с, = (1 +='Н ‘) - с”. Данная дисперсия связывает неопределенность модели 

с, и неопределенность данных с”. Результатом подстановки наблюдаемого изменения 


давления у, в (7) будет вероятность наблюдения у,., в момент х,., для данной модели 


п+1 
на основе первых и наблюдений. 

На рис. 1 представлены прогнозные распределения вероятностей для двух моделей, 
которые обозначены как «Модель 1» и «Модель 2». Видно, что вероятность у,„., для пер- 
вой модели выше, чем для второй. Поэтому «Модель 1» более адекватна, чем «Модель 2». 


Уравнение (7) указывает на то, что вероятность наблюдения у,., в момент х,., будет тем 


п-+1 


выше, чем меньше общая прогнозная дисперсия {ей + с, И наблюдаемое изменение 


давления у, ближе к прогнозируемому (ожидаемому) значению у,„.|. 
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Плотность 
вероятностей 


У 
а) Одинаковые ожидаемые значения, разные дисперсии 


Плотность 
вероятностей 


Ул 
6) Разные ожидаемые значения, одинаковые дисперсии 


Плотность 
вероятностей 


Ул 
в) Разные ожидаемые значения, разные дисперсии 


Рисунок 2 — Три возможных случая прогнозных распределений 
вероятностей для двух моделей 


® Смоделированные данные ГДИС 


Модель 1 
Е: Модель 2 


Изменение давления и его производная, МПа 


10° 10 10° 10 10? 
Время, ч 


Рисунок 3 — Окончательное сравнение Модели | и Модели 2 


На рис. 2 показаны три возможных случая для двух моделей. Так, на рис. 2(а) 
ожидаемые изменения давления для двух моделей совпадают, а общие прогнозные 
дисперсии различаются. В таком случае вероятность наблюдения у. для модели с 
меньшей дисперсией будет выше. На рис. 2(6) ожидаемые изменения давления для 
двух моделей различны, а общие прогнозные дисперсии одинаковы. В этом случае 
вероятность наблюдения у, будет выше для модели, у которой ожидаемое измене- 
ние давления ближе к его действительному значению. На рис. 2(в) ожидаемые изменения 
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давления и общие прогнозные дисперсии отличаются. Тогда вероятность наблюдения 
у,. будет выше для модели, которая имеет меньшую дисперсию и у которой ожидае- 


мое изменение давления ближе кего действительному значению. 

Однако в тех случаях, когда реальное изменение давления ближе к ожидаемому 
значению для модели с болышей общей прогнозной дисперсией, этой модели будет 
соответствовать более высокая вероятность наблюдения у,., в момент х,,,. Это говорит 


п-+1 
о том, что Данная модель, имеющая на данной стадии большую общую прогнозную 
дисперсию, может оказаться корректной моделью в будущем. То есть последовательная 
процедура используется для того, чтобы проверить, насколько верно это предположение. 


0.1 - ` Модель 1 В = 
\. п жж Модель2 / 


Нормализованная суммарная вероятность 


10" 10 
Время. ч 


Рисунок 4 — Нормализированные суммарные 
вероятности 


Модель 1 


Дисперсия, МПа? 


10 10 107 
Рисунок 5 — Оценка дисперсии о” 


Другими словами, подстановка действительного изменения давления в формулу 
для прогнозного распределения вероятностей для каждой модели является процессом 
принятия решения о том, какая из моделей наиболее адекватна на текущей стадии. 
При этом последовательная процедура представляет собой аккумулирование резуль- 
татов принятия решений на всех стадиях. 


Практические рекомендации по использованию метода 


Выбор того, где должна начинаться процедура, является довольно сложной проб- 
лемой. Необходимо принимать во внимание два следующих аспекта из практики ГДИС. 
В течение режима влияния ствола скважины практически отсутствует фильтрация через 
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вскрытую поверхность забоя и стенок скважины, в результате в пласт не поступает входной 
сигнал. Поэтому в течение этого периода ГДИС по результатам измерений давления 
невозможно получить какую-либо информацию о пласте. С другой стороны, как только 
входной сигнал достиг пласта, выходной сигнал теоретически несет всю информацию о 
пласте: средняя проницаемость, скин-фактор, граничные эффекты, влияние неоднород- 
ности и т.д., поскольку изменения давления распространяются внутри пласта с конечной 
скоростью ввиду их диффузионной природы. Для учета этих двух аспектов начальная точка 
выбиралась сразу же после горба кривой, возникающего из-за влияния ствола скважины, на 
графике производной функции давления. 

Теоретически можно выбирать временные шаги в любой последовательности. 
Однако в данной статье временные шаги выбираются в хронологическом порядке. Гра- 
ничные эффекты проявляются, как правило, через определенное время и в процессе 
изменения внутрипластового давления становятся все более и более заметными. Хроно- 
логический выбор временных шагов делает метод более чувствительным к обнаружению 
этих эффектов. 


Модель 1 
ий Модель 2 


10” 


10 10' 10 
Врема, ч 


Рисунок 6 — Значения 2'Н '2 


10 т 


Модель 1 
пееееннние Модель 2 


ю 


=: 


“--.. 
тонны 


Общая прогнозная дисперсия, МПа? 


10 10' 10 
Время, ч 


Рисунок 7 — Общая прогнозная дисперсия 
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0.1 т 


Оценка проницаемости, мкм? 


0.02 - - 
0.01 Е Модель 1 и 
ны Модель 2 
0 т 1 2. 1 1 Ш ё 1 к. 2. 1 4 Ш 
10° 10' 10? 
Время. ч 


Рисунок 8 — Оценка проницаемости 


Кроме того, такой выбор согласуется с концепциями радиуса исследования и времени 
установления. Данный подход также наглядно показывает, каков характер отклонений 
значений давления для наблюдаемых данных от данных, полученных с помощью модели, 
зависящей от времени. Подобная информация очень полезна для принятия решения о том, 
подходит ли модель или нет. 

Хотя обычно за раз прогнозируется только одна точка наблюдения, теоре- 
тически возможно предсказывать любое количество изменений давления в любой 
момент времени на базе фиксированного числа точек данных. Другими словами, с 
помощью первых ип наблюдений за раз можно прогнозировать изменение давления 


не только в момент х,„.,, но также и в момент х,.,. Однако в работе на каждом шаге 


п+1 2 
процедуры делался прогноз только для одной точки данных. 

На каждом шаге процедуры для прогнозирования изменения давления на следую- 
щем временном шаге использовались все данные, начиная с первой точки и заканчивая 
текущей исследуемой точкой. Каждый раз, когда процедура повторяется, количество 
используемых точек данных увеличивается на единицу. Оценки параметров обновляются 
на каждом временном шаге. Новые данные о давлении содержат дополнительную инфор- 
мацию о параметрах при условии, что модель корректна. В качестве алгоритма нелинейной 
регрессии для оценки параметров применялся метод Гаусса — Маквардта со штрафными 
функциями и специальной техникой интерполяции и экстраполяции. 

В начальной точке нет никакой информации о том, какая из моделей лучше. 


Если число возможных моделей обозначено как и, , вероятность, связанная с каждой 


ш 2 


моделью, устанавливается равной |/и.,. 


Для проверки адекватности модели на всем временном интервале должна 
использоваться суммарная вероятность. Суммарная вероятность для каждой модели 
может быть легко вычислена путем последовательного умножения условных вероят- 
ностей на каждом шаге прогноза. 


«Штучний 1нтелект» 172009 321 


Сидельников К.А., Денисов С.В. 


6С 


Из определения условной вероятности следует, что 


ИИ ОВ СЕ (8) 
Ртов( А) 
Поэтому 
_ РО авы У,» У 
РОВ (на у,)= ОЕ (9) 
Это приводит к тому, что 
РОВ а и 5 РОВ а У у,)= 
В РтоБ(у!,..., и а , РгоБ (у... во У) и (10) 
РОВ ны РОВ ыы 
И, 
в РтоБ(у,».... у,) | 


Доверительный интервал, % 


| Модель 1 
ки Модель 2 


10" - 
10 10 10 
Время, ч 


Рисунок 9 — Доверительный интервал оценки проницаемости 


То есть произведение условных вероятностей представляет собой суммарную веро- 
ятность наступления событий у,..., У,, Уи» У,.> В рамках данной модели. Суммарная 


вероятность характеризует степень адекватности модели на всем временном интервале. 

Суммарные вероятности, связанные с каждой моделью, нормализуются так, что 
в сумме они дают единицу. Когда в течение последовательной процедуры нормали- 
зованная суммарная вероятность модели стремится к единице, то можно сделать вывод 
о том, что эта модель является наиболее подходящей. В случаях же, когда модель не 
имеет преимуществ по сравнению с другими моделями, то возможны два варианта: 
следует использовать другую модель или исходных данных недостаточно, чтобы 
провести различие между моделями. 

После п наблюдений суммарные вероятности, связанные с и, моделями, есть 


п, П.,..., П"_. Суммарные вероятности, вычисленные на основе п+1 наблюдения, 


п, 
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есть П!“", П",..., П»'. Нормализованная суммарная вероятность, соответствующая 
7-й модели, определяется как 


П”Р”"*" 
п” —: = =: (1 1) 


У Пт 2 
п рп 
Упр 
= 
где ео — прогнозная вероятность и+1-го наблюдения для ]-й модели, вычислен- 


ной на базе И наблюдений: 
РАН = РОБ (у у, }. (12) 


Оценка скин-фактора, отн. ед. 


Э | Модель 1 


аа Модель 2 
1 4. 


10' 10 
Время. ч 


Рисунок 10 — Оценка скин-фактора 


Выбросы могут привести к чересчур малой вероятности, поскольку прогнозное 
давление будет находиться далеко от резко выделяющегося значения эксперименталь- 
ной величины 8. Однако трудно определить, какие точки данных являются выбросами 
до того, как будет выбрана истинная модель пласта. В результате в работе устанавли- 
вался нижний предел вероятности, что автоматически исключало из рассмотрения 
менее вероятные точки данных. На практике в случаях, когда вычисленное значение 
прогнозной вероятности лежит вне 95 % интервала, оно заменяется нижним пределом, в 
качестве которого выбирается величина вероятности на границе этого интервала. Мате- 
матически это выражается следующим образом: 


РБ, —96. с’ +0, 


5 еслиу, 1 < —1,96. О. +0}; 


Роу, = Роу, +196: о’ +07 


5 еслиу, >, +Ъ96. 0'+0#; (13) 


Роу, и, если, —1,96. 6+0} <у, 151+: 07 +0%. 


Эта модификация делает метод менее чувствительным к выбросам в данных. 

Когда параметры вычисляются на очередном шаге по времени, оценочные 
значения параметров на предыдущем шаге используются как начальное приб- 
лижение. Другими словами, когда параметры определяются на основе п+1 точки 
данных, оценочные значения параметров, рассчитанные с помощью и точек данных, 
выступают в роли начальных оценок. 


«Штучний 1нтелект» 172009 323 


Сидельников К.А., Денисов С.В. 


6С 


Пример дискриминантного анализа двух моделей 
пласта-коллектора 


Модель 1 
анаы Модель 2 


Доверительный интервал, отн. ед. 


10 
Время, ч 


Рисунок 11 — Доверительный интервал оценки скин- 
фактора 

Данные об изменении давления были получены на основе модели пласта с 
непроницаемой внешней границей, к которым добавлялись случайные ошибки, 
распределенные по нормальному закону с дисперсией 2,5.10° МПа’. Информация 
о пласте и насыщающем его флюиде взята из табл. 1. Истинные значения параметров 
К =0,05 мкм”, 5=10, С=0,2 м’ /МПа и г = 600 м. Общее количество точек дан- 
ных составило 81 точку. 


Таблица 1 — Информация о пласте и насыщающем его флюиде 


Радиус ствола скважины (7. ) м 0,1 
Мощность пласта (й) м — 
Объемный коэффициент ( В. ) [1%] / м 1 
Вязкость (и) Па-с 103 
Пористость (ф) 0,2 
Первоначальное давление (Р) МПа 20 
Общая сжимаемость ( с, } МПа ' 10“ 
Рабочий дебит (4) м? /сут 100 
Таблица 2 — Конечные оценки параметров 
Параметры а Модель 1 Модель 2 а 
| мкм” 2,7094 107 4,9720 10? 5.107 
К) 0,2105 .10' 0,9901.10' 1.101 
С м/МПа 1,8734 .107 2,0000 .10`"' 2.10" 
и. м - 6,0050 10? 6.10? 


Для простоты были выбраны две модели-кандидаты: модель радиального 
течения в бесконечном пласте и модель пласта с непроницаемой внешней границей, 
которые далее обозначены как «Модель 1» и «Модель 2» соответственно. 
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Рисунок 12 — Оценка коэффициента накопления 


Конечные оценки параметров приведены в табл. 2. Подобранные по этим 
параметрам окончательные кривые изменения давления и его производной изобра- 
жены на рис. 3. 

Как следует из рис. 4, в течение времени примерно от 0,3 до 20 часов 
нормализованная суммарная вероятность для Модели | практически равна единице, 
т.е. Модель 1 заметно лучше соответствует данным на этом участке по сравнению с 
Моделью 2. Однако после примерно 20 часов нормализованная суммарная вероят- 
ность для Модели 2 быстро возрастает до единицы и остается такой вплоть до конца 
ГДИС. Другими словами, для адекватного описания пласта в его модели необходимо 
учитывать наличие непроницаемой внешней границы. 


ло" 


Доверительный интервал, % 


Модель 1 
ыы Модель 2 


10 10" о 
Рруьаяа. а 


Рисунок 13 — Доверительный интервал оценки коэффициента накопления 


С физической точки зрения, до того как проявятся граничные эффекты, дан- 
ных о наличии внешней границы недостаточно, поэтому Модель 1 более адекватно 
отражает изменение давления. В результате, МИПВ сначала отдает предпочтение 
модели течения в бесконечном пласте, а затем, когда информации становится 
больше, — модели пласта с непроницаемой внешней границей. Такой подход вполне 
согласуется с концепцией радиуса исследования. 

Результаты данного эксперимента наглядно демонстрируют одно из требова- 
ний МПТ, заключающееся в том, что процедура дискриминантного анализа моделей- 
кандидатов пласта-коллектора должна повторяться до тех пор, пока не закончатся 
данные. Если закончить исследование раньше, то возможна ситуация, когда Модели 
1 все еще будет соответствовать значение нормализованной суммарной вероятности 
около единицы, тогда как очевидно, что эта модель не в состоянии объяснить резкий 
рост производной функции давления в виде прямой линии с единичным наклоном в 
конце ГДИС. Однако следует заметить, что нужно анализировать не только конеч- 
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ные значения нормализованных суммарных вероятностей, но также и их динамику 
(например, наличие трендов в течение определенного времени) на всем временном 
интервале ГДИС. 

1200 |- - 


900 |- И: | - 


600 |--- еее еее нениння=# я еее еее ееннннн- + 


Оценка расстояния, м 


300 | - 


№ 


Рисунок 14 — Оценка расстояния до границы 


На рис. 5 оценка дисперсии с” для обеих моделей в начальные моменты времени 
имеет значение, близкое к истинному 2,5.10° МПа’, а после примерно 12 часов 
дисперсия у Модели 1 быстро увеличивается, а у Модели 2 остается примерно на том 
же уровне. 

Величина #'Н\ для Модели 2 всегда выше, чем для Модели 1, что видно на 


рис. 6. Это явилось результатом добавления еще одного параметра (и. ) в модель пласта. 


Общая прогнозная дисперсия (1 + "Н"®) -о” для Модели 1 сначала значительно 


меньше, чем для Модели 2 (рис. 7). Однако с ростом влияния внешней границы наличие 
дополнительного параметра 7’ у Модели 2 улучшает ее согласие с данными, что 


перевешивает вклад этого параметра в увеличение величины #'Н'®. В результате общая 
прогнозная дисперсия для Модели 2 становится меньше, чем для Модели 1. 
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Рисунок 15 — Доверительный интервал оценки расстояния до границы 


На рис. 8 представлены оценки проницаемости, полученные по обеим моде- 
лям. Флуктуации значений оценок в первые моменты времени вызваны нехваткой 


данных. Затем оценки принимают значения, близкие к истинному 0,05 мкм? ‚ а после 


примерно 12 часов оценка проницаемости Модели 1 начинает сильно отклоняться от 
истинного значения из-за возникновения граничных эффектов. Соответствующие 
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доверительные интервалы (рис. 9) сначала уменьшаются с увеличением количества 
данных (причем у Модели 1 доверительные интервалы уже, чем у Модели 2), а 
затем, когда Модель 1 уже не согласуется с данными, ее доверительные интервалы 
начинают резко возрастать на фоне продолжающегося монотонного убывания значе- 
ний доверительных интервалов Модели 2. 

Оценки скин-фактора и соответствующие доверительные интервалы Модели | и 
Модели 2 на рис. 10 и 11 демонстрируют схожее поведение, что и для проницаемости. 
Это связано с тем, что в общем случае между проницаемостью и скин-фактором наблюда- 
ется сильная корреляция. 

Коэффициент накопления, вычисляемый по Модели 1, имеет интересную особен- 
ность (рис. 12 и 13). Теоретически после окончания влияния ствола скважины данные 
ГДИС не содержат больше никакой информации об этом коэффициенте. Поэтому по идее 
оценка коэффициента накопления должна оставаться примерно постоянной. Однако в 
результате сильных граничных эффектов Модель | становится неадекватной наблю- 
даемым изменениям давления, поэтому алгоритм нелинейной регрессии при попытке 
подгонки кривой давления пытается «скомпенсировать» увеличение дисперсии за счет 
изменения всех параметров, в том числе и коэффициента накопления. Все это приводит к 
тому, что и оценки коэффициента накопления, и соответствующие доверительные интер- 
валы, полученные по Модели | на момент окончания ГДИС, трудно назвать приемлемыми. 

На рис. 14 представлены оценки расстояния до внешней границы. До 12 часов 
их отклонение от истинного значения достигает почти 100 %. Однако в конечном 
итоге с увеличением количества данных оценки расстояния сходятся к реальной 
величине 600 м. Аналогично соответствующие доверительные интервалы на рис. 15 
в начальные моменты времени имеют чрезвычайно большое значение (порядка 


10° -+10*). В дальнейшем с ростом влияния внешней границы доверительные интер- 
валы сужаются до приемлемых значений. Если сопоставить рис. 9 и 15, то можно 
отметить интересную особенность. Как известно, возникновение граничных эффек- 
тов приводит к тому, что в пласте устанавливается псевдо-стационарное состояние 
фильтрации флюидов. В результате на изменение давления сильнее начинают влиять 
размеры пласта, а не его проницаемость. Но поскольку при псевдо-стационарном 
режиме фильтрации новой информации о проницаемости меньше, то скорость спада 
кривой доверительных интервалов для проницаемости в Модели 2 также снижается. 
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Метод дискримнантного анал!зу моделей пласта-колектора за результатами пдродинам!чних 
дослджень свердловин 

У стати описуеться новий юльюсний метод для дискрим1нантного анал!зу моделей, який був названий 
методом посллдовного прогнозування Имоврностей. Цей метод базуеться на байесвському висновку. 
Метод е безпосередн!м продовженням способв, що використовують концепщю дов1рних 1нтервалв, але в 
той же час позбавлений 1х основних недолиюв. 
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